Nonlinear delocalization on disordered Stark ladder 
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Abstract. - We study effects of weak nonlineary on localization of waves in disordered Stark 
ladder corresponding to propagation in presence of disorder and a static field. Our numerical 
results show that nonlinearity leads to delocalization with subdiffusive spreading along the ladder. 
The exponent of spreading remains close to its value in absence of the static field. 



Anderson localization leads to suppression of diffusive 

. propagation of linear waves in systems with disorder [1] . 

' In one and two dimensions all states are exponentially lo- 
calized [2,3]. For classical waves nonlinearity is naturally 
present and it is important to understand how it affects 
localization in a random media [4] . At first glance it seems 
that a spreading in space leads to an effective decrease of 
nonlinearity and hence persistence of localization [5] . On 
the other hand it was argued that nonlinear resonances 

, remain overlapped and localization is destroyed by a mod- 
erate nonlinearity which leads to a subdiffusive spreading 

, in space at asymptotically large times [6] . 

Recent experimental progress with nonlinear photonic 
' lattices [7,8] and Bose-Einstein condensates (BECs) in op- 
tical lattices [9, 10] with disorder generated a renewal of 
significant theoretical interest to this problem (see [11-18] 
and Refs. therein). Similar type of problems appear 
also for energy propagation in disordered molecular chains 
[13, 19] that enlarge a field of possible applications. In ad- 
dition the problem of interplay of localization and nonlin- 
earity represents an interesting mathematical problem of 
stability of pure point spectrum with respect to nonlin- 
ear perturbations which led to recent mathematical stud- 
ies [20,21]. 

The numerical studies are mainly done for the discrete 
Anderson nonlinear Schrodinger equation (DANSE) show- 
ing that the wave packet width An spreads at large times 



t in a subdiffusive way with (Art)^ oc t°- and an expo- 
nent a « 0.3 — 0.4 for system dimension d — 1 [6, 14, 17] 
and a « 0.25 for d = 2 [18]. The theoretical estimates 
give a = 2/5 [6] and a = 1/4 [18] respectively. A notice- 
able difference between the theory estimates and numeri- 
cal value a « 0.3 for d = 1 is argued to be related with a 
specific properties of Id Anderson model [18] but further 
clarifications of this point are required (see e.g. [17]). 

In this work we address a new type of question for the 
DANSE model: how a static field force affects the prop- 
erties of nonlinear delocalization? Such a force is experi- 
enced by BECs in a gravitational field or effectively in a 
magnetic field gradient. It can be also effectively created 
by an acceleration of the optical lattice as a whole. This 
creates an effective Stark ladder which already has been 
realized in experiments with cold atoms [22]. In absence 
of nonlinearity a weak static field does not significantly af- 
fect the localization while at strong fields the localization 
length is significantly reduced since less states are energet- 
ically available for hopping over the ladder (see e.g. recent 
studies [23] and Refs. therein). It is not so obvious what 
are the effects of nonlinearity in such a system since the 
nonlinear term is local and is small compared to an energy 
variation for large displacements along the ladder. 

To answer the above question we study numerically the 



I.Garcia-Mata, D.L.Shepelyansky 



Stark DANSE model described by the equation 

dip 

„+l + V'n-l) , 
(1) 

where / is the strength of the Stark field. At / = the 
model is reduced to the usual DANSE studied recently 
in [13,14,16,17]. We fix the units as = ?i = 1 and choose 
a typical set of parameters used here as W = 4, /3 = 1. 
On-site energies En are randomly and homogeneously dis- 
tributed in the interval -W/2 < < W/2. Then the 
localization length in the middle of energy band is £ « 6 
at / = 0. The numerical integration was done by the split 
operator scheme described in [18]. Such a symplectic in- 
tegration with an integration time step At = 0.1 and 0.01 
gives the energy conservation with accuracy 3% and 1% 
for t < 10^ in a presence of strong field / < 2. The total 
number of states was fixed a,t N = 256, we used averag- 
ing over Nd = 15 disorder realisations. The finite value of 
At generates high frequency equidistant harmonics with 
frequency spacing 27r/At. At / = these frequencies are 
located outside of energy band while at / > they, in prin- 
ciple, may give resonant transitions. However, at small At 
the distance between such resonant states is much larger 
than the localization length £ and the matrix elements in 
such cases are exponentially small and do not affect the 
behavior of the system with variation of At (see Fig.l, in- 
set). The integration scheme conserves exactly the total 
probability. 
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Fig. 1: Dependence of the second moment (An)'^ of probability 
distribution on time t for various values of static field / and one 
fixed disorder realisation. Curves from top to bottom at t = 
10* are for / = 0., 0.1, 0.25, 0.5, 1., 2. a.nd p = 1,W = 4, N = 
256, At = 0.1. Inset shows data for / = 0.5 obtained with 
integration steps At = 0.1 (solid curve) and At = 0.01 (dotted 
curve). To suppress fluctuations time average was made on 
logarithmic scale inside intervals A(logiQt) = 0.1. Initial state 
is one lattice site n = with energy in the middle of the band. 



A typical dependence of the second moment (An)^ of 
probability distribution on time t is shown in Fig. 1 
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Fig. 2: Probability distribution Wn = I'i/'nP over ladder site 
n at time t = 10* for / = 0., 0.1, 0.25, 0.5, 1., 2. (curves from 
small (inside) to large (outside) values of \n\). Parameters are 
the same as in Fig. 1. 



for various values of /. The distribution of on-site prob- 
abilities Wn at a final time t = 10^ is shown in Fig. 2. 
The fits of data of Fig. 1 show an unlimited subdiffusive 
growth (An)2 cx t« with the exponent a = 0.275 ± 0.006, 
0.291 ± 0.005, 0.276 ± 0.006, 0.272 ± 0.004, 0.269 ± 0.016 
for / = 0., 0.1, 0.25, 0.5.1. respectively. The fits are done 
in the interval 5 < logj^Q t < 8 for the last case and 
3 < logj^Q t < 8 for all other cases. The data show no 
significant variation of a with / oven if one realization 
at finite time may have noticeable fluctuations of a being 
larger than a formal statistical error (sec below). At large 
values of / = 1., 2. the localization length i becomes rather 
small (it can be estimated as w y^(An(t = 1000))^) and 
during a long time interval there is not spreading over the 
ladder. For / = 1. the growth appears at t > 10'"'. There 
is no visible growth for all computational times for f — 2. 
We interpret this as very low transition rates over local- 
ized states in the case of small localization length £ 1 
at / = 2. It remains unclear if localization persists or 
disappears for such small £ at very large times. For the 
cases with clear delocalization at / = 0., 0.1, 0.25, 0.5 the 
distribution of w„ over n have a form of homogeneous 
"chapeau" centered near the initial state n = 0; its width 
grows with time, approximately in agreement with the sec- 
ond moment growth (see Fig. 2). 

To obtain more statistics we perform averaging over Nd 
disorder realisations. The data are presented in Fig. 3 
for the second moment and in Fig. 4 for the probability 
distribution for = 0;1. At / = we obtain the ex- 
ponent a = 0.302 which is comparable with the values 
a ~ 0.33 found in previous studies [14, 17, 18]. The value 
of a decreases by about 10% when the static field is in- 
creased up to / = 0.5. This decrease is well visible even 
if it is not very large and is comparable to the statistical 
variations of a at / = discussed above. We also com- 
puted the dependence on time for the participation ratio 
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^ = 1/ < J2n >■ It be characterized by the depen- 
dence^ cx f with the exponent = 0.120±0.001(/ = 0.5), 
0.131±0.001(/ = 0.25), 0.159±0.002(/ = 0). These values 
are compatible with the usual relation ly = a/2. 



since they are separated by many localization lengths £ of 
the linear problem. Nevertheless the propagation on these 
far ends goes in a correlated way since the total energy 
Etot ~ J2n /'^'"^n is exactly conserved. 




Fig. 3: (Color online) Dependence of the second moment (An)^ 
on time t, average is done over Nd = 15 disorder realisa- 
tions. The curves from top to bottom at t = 10* are for 
/ = 0., 0.25, 0.5 at /3 = 1 and / = 0.5 at /3 = 0. The fits, 
shown by thin straight lines, give the exponent of growth at 
(3=1: a = 0.302 ± 0.001(/ = 0) , 0.262 ± 0.001(/ = 0.25), 
0.241 ± 0.002(/ = 0.5). Inset shows the dependence of par- 
ticipation ratio ^ on time for /3 = 1,/ = 0.5 (thick curve), 
the straight line shows the fit dependence with the exponent 
f = 0.120 ± 0.001. Other parameters are as in Fig. 1. 

The averaged probability distributions for the cases of 
Fig. 3 at / = 0.5 are shown in Fig. 4 at different mo- 
ments of time. For (3 — the distribution is localized 
being frozen in time. It drops faster than exponential due 
to the presence of static field showing a qualitative dif- 
ference between Stark localization and usual exponential 
Anderson localization. For (3 = 1 the probability spreads 
over the whole lattice forming a homogeneous plateau in 
the center. The interesting property of this distribution 
is its approximate symmetry with respect to the initial 
state n = 0. It is clear that the conservation of energy 
imposes such a symmetric spreading. Indeed, the width 
of the plateau is approximately Sn w 80 and the energy 
variation on such a distance is SE « fSn « 40 that is 
much larger than the energy band B Ki 6 at f = 0. Due 
to energy conservation at |/| > the spreading can con- 
tinue unlimitedly only in approximately symmetric way. 
This excludes the possibility of a compact packet which 
moves over a lattice on larger and larger distances in some 
stochastic way (as discussed in [13]). A quasi-symmetric 
spreading seen in Fig. 4 looks rather natural in view of 
total energy conservation at |/| > 0. However, it raises 
an interesting problem of statistical entanglement of prob- 
abilities Wn on opposite ends of the plateau. Indeed, on 
such a distance the probabilities seems to be uncorrelated 
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Fig. 4: (Color online) Averaged probability distribution Wn as 
a function of the lattice site n for the cases of Fig. 3 at / = 0.5 
(average is done over the same 15 disorder realisations). Top 
panel is for /3 = and bottom panel is for /3 = 1. Curves are 
drown from inside (small |n|) to outside (large |n|) for t = 10^ 
(magenta), 10* (blue), lO** (red) and 10* (black). For (3 = 
(top panel) the state is almost the same for all times from 
t = 10^ tot = 10*. 



In [6, 14, 18] it was argued that an infinite spreading is 
possible since the nonlinear frequency shift Suj ~ P/An re- 
mains comparable with the frequency spacing Auj ~ 1/ An 
between frequencies of excited An modes. On the Stark 
ladder this condition seems to be violated at large An 
since Aw ~ / 3> 5w ~ /3/An. However, the situation is 
more subtle. Indeed, we can write DANSE (1) in the basis 
of linear eigenmodes. The time evolution amplitudes Cm 
in this basis is described by equation (see e.g. [6, 18]): 

g(j ^ 

(2) 

where m marks the center of eigenmode inside the ladder 
and eigenenergies are randomly distributed inside the 
energy band width of approximately the same size _B ~ 4 
as at / = 0. Vmm' mim[ l^^^"^ are the transition matrix 
elements induced by nonlinearity. From this equation it is 
clear that 4-waves resonance conditions are satisfied if the 
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frequency detuning Aa;4 of these 4-waves remains small: 
Aa;4 = f{m + mi-m' -m[)+em + Cmi - em' - £m[ < 

(3) 

Thus the spreading over the ladder can proceed only over 
such modes where m + mi — m' — m[ = and thus 
Aoji = Cm + Cmi — Cm' — ^m[- Siuce all Cm are inside 
the frequency band B it is possible to have Alo/^ ^ 1/An 
so that the resonant detunings will remain small compared 
to nonlinear shift Sui ~ /3/An even at large An values. 

The rate of spreading is still determined by the same 
estimates as in [6, 14, 18] since we still have dC/dt ~ (3C^ 
and the theoretical exponent is a = 2/5 being independent 
of / if /£ < 1 so that the local transition rates remain the 
same as at / = 0. The numerical results presented here 
show weak dependence of a on / for f < l/£ being in a 
satisfactory agreement with this theoretical estimate. The 
deviation of a « 0.3 from the theoretical value 2/5 still 
should be better clarified both for / = and |/| > 0. 

In summary, we demonstrated that, in presence of a 
static field applied to a lattice with disorder, a nonlin- 
earity still produces complete delocalization with a sub- 
diffusive spreading over the ladder. The exponent of the 
spreading remains close to the value without the force. 
The spreading forms a homogeneous distribution of prob- 
ability inside a certain plateau. Due to the conservation 
of total probability and energy the far away parts of this 
plateau remain statistically entangled even being a large 
distance apart from each other. The obtained results can 
be tested in experiments with nonlinear photonic lattices 
and BEG atoms in optical lattices with a static field. 
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